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Current efforts in observational cosmology are focused on characterizing the 
mass-energy content of the Universe. We present results from a geometric 
test based on strong lensing in galaxy clusters. Based on Hubble Space Tele- 
scope images and extensive ground-based spectroscopic follow-up of the mas- 
sive galaxy cluster Abell 1689, we used a parametric model to simultaneously 
constrain the cluster mass distribution and dark energy equation of state. 
Combining our cosmological constraints with those from X-ray clusters and 
the Wilkinson Microwave Anisotropy Probe 5-year data gives = 0.25 ± 0.05 
and = —0.97 ± 0.07 which are consistent with results from other methods. 
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Inclusion of our method with all other techniques available brings down the 
current 2a contours on the dark energy equation of state parameter Wy, by 
about 30%. 



Measurements of the Hubble diagram for Type la supemovae ([7] |2] |J] |?]) combined with 
constraints from the Wilkinson Microwave Anisotropy Probe (WMAP5) (|5|[6|), cosmic shear 
observations ([^[S][^[7^[77]), cluster baryon fractions ([72]), cluster abundances (13) and baryon 



acoustic oscillations (BAO) from galaxy surveys ([74][7^[76|) suggests that ^ 72% of the total 
energy density of the Universe is in the form of an unknown constituent with negative pressure 
- the so-called dark energy, that powers the measured accelerating expansion. These observa- 
tions probe the equation-of-state parameter w^, defined as the ratio of pressure to energy density, 
through its effect on the expansion history and structure of the Universe. The current goal of 
cosmology is to understand the properties of dark energy by placing tighter constraints on its 
equation of state. In the currently favored flat ACDM model dark energy is attributed to a 
cosmological constant, for which = —1. Type la supernovae, baryon acoustic oscillations, 
cluster abundances and cosmic shear appear to be very promising techniques to tighten con- 
straints on the equation-of-state parameter in the near future. Because all of these techniques 
have biases, systematics and degeneracies, it is only in combination that robust estimates of 
cosmological parameters can be obtained. 

The most recent census from a combination of techniques suggests that ^ 72% of the en- 
ergy density in the Universe is in the form of dark energy that is powering the accelerating 
expansion of the Universe. In the progression towards In this work, we present results from a 
technique that exploits the strong gravitational lensing of distant background galaxies by mas- 
sive galaxy clusters. Through their effect on the local space-time geometry, massive foreground 



^The currently favored 'concordance' cosmological model that best describes the Universe is the ACDM 
paradigm in which the bulk of the matter and energy density are dominated by dark matter and dark energy with 
baryons contributing only ~ 5%. 



structures cause the deflection and shearing of hght rays originating from distant sources. In 
the case of strong lensing, the hght beams are deflected so strongly that they can often result in 
the observation of several distorted images of a given single background galaxy. The positions 
of these multiple images depend strongly on the detailed properties of the lens mass distribu- 
tion ([7^[7S][79|). Because the image positions also depend on the angular diameter distance 



ratios between the lens, source and observer, they encapsulate information about the underly- 
ing cosmology. We capitalize on this dependence on the geometry to derive constraints on the 
cosmological parameters (the mean matter density) and w^. 

Constraining the energy content of the Universe using multiple sets of arcs in cluster lenses 



has been explored in the past (20 21 22 23 24 25 26). In particular, simultaneous inversion of 



the lens and derivation of cosmological constraints can be performed based on the cosmological 



sensitivity of the angular size-redshift relation with sources at distinct redshifts ( [27] ). In this 
method, the angular diameter distance ratios for 2 images from different sources defines the 
'family ratio' S, from the cosmological dependence of which constraints on and are 
extracted: 

where z\ is the lens redshift, z^i and Zs2 are the two source redshifts, and D{zi, Z2) is the angular 
diameter distance. 

Application of this method to the cluster Abell 2218 using 4 multiple image systems at 
distinct redshifts, resulted in < 0.37 and < —0.80 for a flat Universe ( |25| ). A recent 
feasibility study demonstrates that the degeneracies of this technique are entirely distinct from 
those of other cluster methods and that combining the results from several simulated clusters 



with > 10 multiple image families in each can provide a powerful probe of dark energy ( [26| ). 

We have applied this technique to the massive, lensing cluster Abell 1689 at redshift z = 
0.184. Based on images from the Advanced Camera for Surveys (ACS) aboard the Hubble 



space Telescope (HST) this cluster has 114 multiple images from 34 unique background galax- 
ies, 24 of which have secure spectroscopic redshifts (ranging from z ^ 1 to z ^ 5) obtained 



with the Very Large Telescope (VLT) and Keck Telescope spectrographs (|2^|2S|). Abell 1689 
is amongst the richest clusters in terms of the number density of galaxies in its core. It is also 
amongst the most luminous of galaxy clusters in X-ray wavelengths, with an absolute X-ray 



luminosity of = 20.74 x 10^^ W ( |29| ). Observationally, Abell 1689 consists of two groups 
of galaxies : a dominant one located at the center coincident with the peak of X-ray emission, 
and a secondary concentration about 1 arcminute North-East of the main one. Studies have 
shown that this second northern group is at a slightly higher redshift, suggesting thus that these 
two groups might actually be merging The projected mass enhancement produced by such 
a merging configuration, could therefore explain the stunningly large number of multiple im- 
ages identified in this cluster. Previous work shows that the mass distribution of Abell 1689 is 
well modeled with a set of parameterized elliptical pseudo-isothermal lensing potentials (28). 
We utilized the most recent parametric model of Abell 1689, which is able to reproduce the 
observed image configurations to within an average positional accuracy of 2.87 arcseconds, as- 
suming a ACDM cosmology. We solved the lens equation in the source plane for Abell 1689 as 
it is computationally efficient. Inverting in the lens (image) plane provides additional informa- 
tion but is computationally prohibitive at present ( |2g| ). Our simplified model that has a total of 
21 free parameters consists of two large-scale potentials, a galaxy-scale potential for the central 
brightest cluster galaxy (BCG), and includes the modeling of 58 of the brightest cluster galax- 
ies. Therefore, we explicitly include the effect of substructure in the lens plane and assigned 
potentials associated with bright cluster galaxies. The velocity dispersion and scale radii of all 
but one (the BCG) of the cluster galaxies were assumed to follow empirically motivated scaling 
relations, which have been previously utilized to model cluster lenses ( [37] [32] ). 

Despite the large number of multiple images observed, not all of them can be utilized to 
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constrain cosmology. From the initial 1 14 images, we only used those (i) with robust, measured 
spectroscopic redshifts; and (ii) excluded those in the regions of the cluster with low S/N in the 
mass reconstruction. This selection results in the culling of multiple images that lie in the most 
uncertain regions of the mass distribution. Moreover, we identified several bright spots in some 
well resolved multiple images, which we used to increase the number of families. Applying 
criterion (i), resulted in a catalog of 102 images. Imposing criteria (i) and (ii), we obtained a 
catalog of 28 images arising from 12 families all with measured spectroscopic redshifts (Fig.|4]), 
providing a total of 32 constraints. We assume flatness as a prior {Qtot = 1) and fix the Hubble 
parameter at Hq = 74km/s/Mpc as our cosmography test is not sensitive to the value of 
the Hubble parameter. 

For each of the observed image systems with n images, we determine the goodness of fit for 
a particular set of model parameters using a source plane x^, 

1=1 

where is the source plane position corresponding to image i, ^^5^ is the family bary center, 
M is the magnification tensor, and is the total (observational and modeling) error. The total 

was obtained by summing over families and was used in conjunction with a Markov Chain 
Monte Carlo (MCMC) sampler to probe the posterior probability density function (PDF) as a 
function of all relevant model parameters (SOM). The key degeneracies with cosmological 
parameters for this technique arise from the velocity dispersions, ellipticity and core radii of the 
large scale mass clumps in the model (Fig. S4). 

The angular resolution of HST images is on the order of 0.1 arcseconds. However, the 
modeling errors are generally larger due to the complexity of the cluster mass distribution, as 
well as the effect of intervening structures along the line of sight. We quantified the errors 
due to the presence of structure along the line of sight using the Millennium Simulation halo 
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catalogs ( [35| ). We quantified the errors on an image by image basis. By randomly slicing 
through snapshots of the simulation and tiling them at redshift planes between the observer and 
the source, we constructed 1000 line-of-sight realizations. We then ray traced through each 
realization with the Abell 1689 model included to estimate the effect of intervening halos on 
image positions. In most cases, these line-of-sight halos perturbed image positions but did not 
alter the multiplicity of the images. These perturbations induce positional displacements of 
the order of 1 arcsecond. Therefore, about 1 arcsecond of the error between observed image 
positions and model image predictions can be attributed to the presence of structure along the 
line of sight behind the cluster. The presence of projected correlated and associated large-scale 
structure (filaments) increases the cross-section to strong gravitational lensing making these 
clusters more efficient lenses, however, simulations show that this is a sub-dominant effect to 
that of unassociated distant large-scale structure. Therefore, the presence of aligned correlated 
large scale structure at the redshift of the cluster does not scupper the recovery of cosmological 
parameters [5 

Earlier work on reconstructing the mass distribution of several massive lensing clusters, 
shows that the association of dark matter substructures with the locations of the brightest cluster 



galaxies is well matched by that derived from the Millenium Simulation (|5^|J2]). This supports 
our use of luminosity-mass scaling relations to map substructure in the cluster. 

The second potential source of error arises from modeling uncertainties for the substructure 
in the lens plane. This is likely due to scatter in the assumed scaling relations for the velocity 
dispersions and scale radii of cluster galaxies. Although the mean correlations between these 
variables may be well described by simple scaling relations, individual galaxies can deviate 
substantially from them, introducing errors into the parameter recovery. In order to quantify 



^Filaments aligned at finite incHnations behind the cluster (thus breaking azimuthal symmetry) do yield larger 
deviations between the projected and multi-plane models compared to the symmetric case. For the non- symmetric 
case, the deviations are typically on the same level as the errors due to scatter in the cluster galaxy properties. 
These deviations are still subdominant to the effects of uncorrected line-of-sight (LOS) halos. 



these modeling errors, we performed Monte Carlo simulations of the lens system assuming a 
20% scatter in the galaxy scaling relations. This scatter induced modeling errors on the back- 
ground galaxy image positions that in some cases were as large as ^ 1 arcsecond. Therefore, 
the estimated errors from substructure effects in the lens plane and along the line of sight are 
comparable for the selected multiple image families in Abell 1689. 

We use the catalog of 102 images (including images with photometric redshift estimates) 
and our estimates of the observational and modeling errors to obtain the marginalised PDF in 
the — Wy, plane. Adding in quadrature the systematic errors identified above to the posi- 
tional uncertainties, results in a best model for Abell 1689 with a well defined though broad 
degeneracy between and Wy,. The "concordance" model of Vt^n ^0.3 and —1 lies 

within the la contour, but the way the degeneracy is pushed against the prior limits suggests 
either a bias in the mass modeling, or misidentified images (Fig. S5). Our simulations show 
that the errors in photometric redshift determination methods at present are too large and limit 
the efficiency of our technique and introduce biases (see SOM). Therefore, we excluded all 
images with photometric redshifts in our modeling, bringing down the number of image fam- 
ilies used from 34 to 24. Even with this cut, some images were badly reproduced with very 
high RMS positional deviations, in particular, those that lie in complex crowded regions of the 
cluster (regions where the signal-to-noise of the mass map is low due to the presence of several 
bright cluster galaxies in close proximity). These outlier images from 16 families highlight 
complex regions in the lens mass distribution, not handled adequately by our simple parametric 
model. The outlier image systems also tend to have higher redshifts; thus they are likely to be 
more affected by uncertainties in modeling the intervening line of sight structure as well. We 
thus deliberately disregarded images that lie in the most complex regions in the mass distribu- 
tion. While the recovery of cosmological parameters is insensitive to the choice of profile, both 
the PIEMD (Pseudo-Isothermal Elliptical Mass Distribution) and NEW (Navarro-Erenk- White) 
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provide comparable constraints on and Q^n (Fig S3), observationally some clusters are better 
fit to one or the other model. Abell 1689 is best fit with a PIEMD profile. The final culled image 
catalog thus contains 28 images, derived from 12 families at redshifts ranging from zs = 1.15 
to Zs = 4.86, all of which are spectroscopically measured (Table SI and Fig.|4]). 

Optimising our model with all the spectroscopically selected images, including outlier im- 
ages does not result in significant constraints on either and and yields an averaged re- 
duced = 0.08. This indicates that the line of sight and scaling relation errors are likely over- 
estimated. Thus, we optimized again but this time excluding the outlier images. In this iteration, 
we obtained an averaged reduced ^ 28, indicating that now the errors were somewhat under- 
estimated (Fig. [8]). Owing to the large estimates for the modeling errors, constraints obtained in 
this fashion from a single cluster lens are fairly modest. However, combining our results with 
those from X-ray clusters and WMAP5 leads io Q^n = 0.25 ± 0.05 and = -0.97 ± 0.07 
(Fig. [8]), which is consistent with the values derived from combining WMAP5 with SN and 



BAO, = 0.265 ±0.16 ±0.025 and = -0.96 ± 0.06 ± 0.12 (jJ^). 

Our results when combined with the results from WMAP5 the supernovae "Gold 
sample" 0, SNLS project (39) and SNEssence, and the BAO peak from SDSS (|76]), give 
0.23 < < 0.33 and -1.12 < < -0.82 at the 99% confidence level (Fig. [3]). This 
combination of all current viable probes brings down the overall error in by about 30%. 
Therefore, the combination of cluster methods with WMAP5 has comparable potency to the 
combination of other cosmological probes. 
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Figure 1: The critical lines for a source at z = 3 are overplotted in yellow on the HST ACS 
image of Abell 1689. The lensing mass model used is the one from which we derived cos- 
mological constraints. In addition to 2 large scale clumps and the BCG, this model includes 
the contribution of 58 cluster galaxies. The positions of cluster galaxies are marked with green 
crosses. Overplotted in white are the 28 multiple images arising from 12 families that we used 
in this work; the red circles mark the positions of the rejected images. 
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Figure 2: Left) The results from the simultaneous Bayesian optimization of the detailed mass 
distribution and cosmological parameters in the — plane for Abell 1689 using the 28 
multiple images belonging to 12 families at distinct redshifts as constraints from strong tensing 
including only observational errors. The plotted contours are the 1, 2 and 3-a confidence levels. 
(Right) The results from combining cosmological constraints from WMAP5-\- evolution of X-ray 
clusters +cluster strong tensing (cluster only methods); the 1 and 2a contours are plotted, blue 
contours - constraints from WMAP5, pink contours - X-ray clusters, orange contours - cluster 
strong tensing. We multiplied the likelihoods for the various techniques to obtain this plot. The 
degeneracy directions for X-ray clusters and cluster strong tensing are orthogonal. 
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Figure 3: Combination of constraints from strong lensing, the WMAP5 data (IF), the supernovae 
"Gold sample" (4), SNLS project (39) and SNEssence, and the BAG peak from SDSS(16). Con- 
tours for non-SL cosmological probes come from the WMAP plotter. The overplotted contours 
are all 1, 2 and 3-a confidence levels. 
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Supporting Online Materials 
Mass reconstruction of Abell 1689 



As a starting point for the cosmography work, we use the strong lensing based mass model 
presented in Limousin et al. 2007. We summarize this model briefly below. The mass distri- 
bution in Abell 1689 was reconstructed in a parametric fashion using detected strong lensing 
features from deep HSTACS observations and extensive ground based spectroscopy. We find 36 
multiply imaged systems yielding 102 images, of which 24 systems have reliably determined 
spectroscopic redshifts. These systems span a redshift range between z = 1.1 and z = 4.9 and 
provide a total of 136 constraints for the analysis. The multiple images used and their relevant 
properties, including location and redshift are shown in Table 1. The resulting mass distribution 
in this cluster is found to be bimodal: the dominant component is a central large scale mass 
clump/halo coincident with the location of the brightest cluster galaxy (BCG) and the center 
of the X-ray emission. The second mass component is a another large scale mass clump/halo 
in the north-east part of the ACS field, it also contributes substantially to the overall mass dis- 
tribution and is required by the data to reproduce several of the multiple image geometries. In 
addition, this model includes the galaxy scale dark matter haloes associated with the all the 
identified cluster members (totalling 267) detected within the ACS field providing the smaller 
scale structure in the dark matter distribution. The critical lines for this mass model are shown 
in Figure |4} To keep the mass model simple and yet include the contribution of small scale 
structure, standard scaling laws with luminosity are assumed for the cluster galaxies. Thus the 
final mass model has 33 free parameters. The large number of observational constraints used 
(136 constraints) makes Abell 1689 the most reliably reconstructed cluster to date, with a mass 
distribution constrained at the per cent level, a level required in order to exploit Abell 1689 
efficiently as a cosmological probe. Note that this is the baseline input mass model that we 
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started with, eventually, to extract cosmological parameters we pared down the model to 
include only 58 of the cluster galaxies and only 28 images all with measured spectroscopic 
redshifts. Details of the selection are presented in the paper. 



Strong lensing modeling 

Although numerical simulations suggest that dark matter halos are well described by the spher- 
ically averaged universal Navarro-Frenk- White (NFW) density profile, in lensing studies we 
routinely prefer the truncated Pseudo Isothermal Elliptical Mass Distribution (PIEMD) poten- 
tials for the extra degree of freedom in allowing ellipticity and due to fact that these models 
have finite mass. The 3 -dimensional PIEMD density distribution is given by: 

o(r) = — (3) 

The PIEMD potential is parameterised by a central density, po, related to the central velocity 
dispersion of the potential, ctq, and two characteristic radii that define the scales on which the 
slope of the density profile changes. The profile is flat in the inner region, then isothermal 
(p r~^) between Vcore and Vcut^ and steeply decreases (p r~^) beyond Tcut- In addition, the 
ellipticity e and the position angle 9 further characterise the shape of the potential. The PIEMD 
potential has been extensively and successfully used to model galaxy clusters ([7^|7^|2S|). 



In the inner regions of clusters that are pertinent to the calculation of strong lensing effects, 
cluster galaxies are mostly early-type galaxies, and have been found to follow the scaling rela- 



tions like the Faber- Jackson relation ( |?7][?2]|?3| ). Utilizing these observed empirical scalings in 



our models, we assume that the masses of cluster galaxies scale with their luminosity L, via the 
following scaling relations: 
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Figure 4: Mass model for Abell 1689 published in Limousin et al. (2007): the critical lines for 
a source at z = 3 are overplotted in yellow on the HSTACS image of Abell 1689. This is done 
for the best-fit lensing mass model that includes 267 cluster galaxies. Overplotted in red are all 
the multiple images including those with photo-metric redshifts. Cluster members selected in 
Abell 1689 using the color-magnitude relation are marked with green crosses. 
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Figure 5: The identification of multiple images in Abell 1689 from the high resolution HST 
ACS data: shown are a few examples wherein the white circle plotted is the positional accuracy 
of each image. The images are labeled as A.B where A denotes the family number and B the 
image number within family A. Note that even within a given family of images of the same 
source, the error varies due to the variation in S/N of the images determined by the background 
that in turn depends on the location of the images. 
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Table 1: Catalog of 28 images used in this work to constrain cosmology. The columns are 
ID, R.A., Decl., z, x (positional uncertainty in CCD row), y (positional uncertainty in CCD 
column), RMS (error between predicted and observed image positions). All images have secure 
spectroscopic redshifts. 
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where is the luminosity of a typical galaxy in the cluster; Vcore is fixed to an arbitrary small 
value of 0.15 kpc;r^^^ and ctq depend on the properties of the cluster and are left as free pa- 



rameters. The small value of rcore is in agreement with the results of from the modeling 
of individual lensing galaxies wherein it appears that early-type galaxies are isothermal in their 
inner regions. According to these definitions, the halo mass is proportional to aQTcut and the 



mass-to-light ratio is constant (31 ) 



Cosmological angular diameter distances are modeled according to the definition from {45) 

as, 

Da = '^ r dz (nUl + zf + (1 + zf^-^^'^nX'" . (5) 

1 + Z Jzi ^ ^ 

In this equation, we assume a flat Universe (choice of prior). In addition, since we assume 
matter and dark energy to be the only significant components in the Universe, the sum of their 
respective parameters Vl^ + = 1- Finally, we assume the Hubble constant to be = 



74 km/s/Mpc [33), As for the two cosmological parameters (f^m, ^x), we assign the following 
uniform priors Vl^ G [0, 1] and G [—2, 0]. 

In optical images, galaxies in Abell 1689 form 2 groups : one predominantly in the centre, 
and the second located about one arcminute to the North-East. As in ( |2g| ), our final model for 
Abell 1689 consists of 2 cluster-scale halos, one galaxy-scale halo to model the BCG at the 
center of the dominant group. On top of that, we add the color-magnitude selected catalog of 



cluster member galaxies from ( [2g| ), but with an additional cut in deflection angle. Only halos 
that produce a deflection larger than 0.07 arcsecond (i.e. about half the observational uncertainty 
on the image positions) are selected. Thus, we obtain a catalog of 58 cluster member galaxies 
brighter than vtlk < 18.11. Each galaxy is assigned a PIEMD potential and the scaling relations 
outlined above are adopted. In total, our final mass model has a total of 21 free parameters, 
with 124 image constraints. The adopted best estimate values for the mass model parameters 
are reported in Table 2. 
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Table 2: Estimated values and errors at la confidence for the parameters of the mass model of 
Abell 1689; " the position of each mass component is given relative to the cluster centre; ^ the el- 
lipticity corresponds to the lenstool input ellipticity for a PIEMD potential e = (a^ — 6^)/ 
where a and h are major and minor axis of the PIEMD mass distribution; the orientation 9 in- 
creases from West to North. 



Before launching any computationally intensive runs with the full set of images, we perform 
a quick check with a scaled down version of our mass model excluding cluster member galaxies. 
For the two smooth cluster-scale potentials and the BCG potential, we optimise the positions, 
ellipticities, orientations, core radii, cut radii, and velocity dispersions. We also let Vt^ and Wy, 
free to vary between the limits stated above. As constraints, we use two systems at redshifts 
3.050 and 2.533. These are both 5-fold multiply imaged systems with most of their images lo- 
cated far from any bright cluster galaxy. In addition, these 2 systems form two Einstein crosses, 
with most of their images appearing at the same distance from the cluster center. Using these 
two systems, the estimated cluster mass obtained is rather insensitive to our choice of the form 
of the density profile. Consequently, our estimated values for Vt^ and Wy, from these systems 
will be independent of our detailed modeling of the lens potential. We obtain a degeneracy 
between Vt^ and Wy, in agreement with our theoretical expectations ( [26| ). Although we obtain 
large error bars on both Vt^ and Wy,, the favored cosmological parameter values {Q.^ = 0.3 and 
Wy, = —1) are well centred in the la contour. For all subsequent runs, we use the 3a confidence 
intervals of each mass model parameter derived from this run as priors. 
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Note that the choice of parameterization of the mass model (NFW profile versus PIEMD) 
does not impact recovery of cosmological parameters. While NFW profiles are found to be 
a good fit to simulated massive galaxy clusters, in reality the NFW is not always the best fit 
to lensing data of observed galaxy clusters. From our extensive experience is modeling the 
mass distributions of lensing clusters combining constraints from both the strong and weak 
lensing regime, we find that for a majority of clusters the PIEMD (Pseudo-Isothermal elliptical 
mass distribution) is a significantly better fit that the NFW. This is most definitely the case for 
Abell 1689. In our earlier work on modeling the mass distribution of this cluster presented in 
Limousin et al. (2008), we investigated a range of models NFW, generalized elliptical NFW, 
PIEMD and power law models. Ranking model fits using Bayesian evidence derived by using 
the MCMC sampler we found that the PIEMD model offered the best fit to the current data. 
Therefore, as a starting point for the cosmography work we started with the best- fit mass model 
parameters which are of course then simultaneously optimised within reasonable ranges (the 
model is not fixed). 

In order to test the role of this choice of parametric model on the recovery of cosmological 
parameters, we have performed detailed simulations with many realizations of both the NFW 
and PIEMD modeled clusters. As we show in Figure 3 the choice of mass model does not bias 
recovery of input parameters. In the simulations, for both mass distributions, it was assumed 
that we have constraints from 20 image families all with measured spectroscopic redshifts. The 
simulated cluster are taken to be at the same redshift z = 0.2, each with 20 families created from 
the same background source catalog: one is PIEMD and one is NFW. There are no subclumps 
included as the issue to settle here is the dependence of the choice of the smooth model. Typical 
observational errors were assumed for both cases. It turns out that the PIEMD contours are 
slightly wider, probably due to the additional parameter needed. 
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Figure 6: Left Panel: The recovery of input cosmological parameters with a simulated cluster 
modeled with an NFW profile utilizing 20 image families. Right Panel: The recovery with the 
same cluster and background image catalog as above but with a cluster modeled with a PIEMD 
profile. 
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Figure 7: Tabulation of the key parameter degeneracies: cosmological parameters are degener- 
ate with the following lens model parameters 24^elocity dispersion, ellipticity and core radius. 
This is shown in the results from the combined Bayesian optimization of the mass distribution 
and cosmology in Abell 1689. In this plot, emass is the ellipticity of the mass; rcore = 
the core radius and sigma the velocity dispersion. The 2 large-scale clumps and the BCG are 
denoted by 1, 2 and 3 respectively. 



Inventory of systematic errors in the modeling of Abell 1689 

Due to the large number of images found in Abell 1689, we can expect the constraints on 
cosmology to be very tight. However, it is known that parametric models are not flexible enough 
to reproduce all the images at the observational level of accuracy (|?^|?^. Moreover, an initial 
analysis with the 0.13 arcsecond positional uncertainty found above, yielded different modes 
in our parameter space with similar likelihoods that provided very different estimates for 
and w^. This often happens when systematic errors are not adequately taken into account. 
Therefore, we present a preliminary analysis of the systematic errors for Abell 1689 in this 
work. We treat two sources in this work: systematic errors arising from cluster galaxies in the 
lens plane and halos along the line of sight. 

Substructure in the lens plane: The velocity dispersions and scale radii of cluster member 
galaxies are likely to display significant scatter about the scaling relations we assumed to 
model them. If not properly accounted for, this scatter can introduce biases into the parameter 
recovery. In order to quantify this effect on an individual image basis, we perform Monte 
Carlo simulations of the lens system. We first obtain a set of source locations by mapping each 
observed image back to the source plane using the initial cluster model. We randomly draw ctq 
and Tcut parameters for the cluster galaxies from Gaussian distributions with mean values 
obtained from the above scaling relations and standard deviations equal to 20 per cent of the 
mean. We then lens the source locations to the image plane using the cluster model with each 
realization of the galaxy population. We find that a 20 per cent scatter in the velocity 
dispersions and scale radii can induce modeling errors on the image positions that can be as 
large as ^ 1 arcsecond. Put in another way, in trying to avoid potential biases in the 
cosmological parameters, we should not expect our simplified models to reproduce image 
positions to within ^ 1 arcsecond given the limitations on modeling the galaxy population. 
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Figure 8: The results from the simultaneous Bayesian optimization of the detailed mass distri- 
bution and cosmological parameters in the Q^n — plane for Abell 1689 including all available 
families of multiple images as constraints including all sources of error from observations, the 
line of sight structure as well as errors in scaling relations employed to model cluster galaxies. 
For this iteration all the 102 images were used including 58 cluster galaxies, error from all sys- 
tematics were added to each image. The RMS error in reconstructing image positions was large 
signaling that the mass model was not accurate enough. 
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Substructure along the line of sight: We also quantify the errors due to halos along the line 
of sight using the Millennium Simulation (MS) halo catalogs ( |35| ), which are freely available 
for download. We construct lens planes by selecting randomly oriented slices from each MS 
snapshot. We obtain a catalog of halo positions and masses (down to 10^^ Mq) in each slice 
and project their positions along one direction. We then place analytic NFW potentials at each 
halo location, using the mass and redshift to assign a concentration parameter through scaling 
relations obtained from N-body simulations in the literature. The lens planes are tiled along 
the line of sight and the cluster model is inserted at the appropriate redshift. We ray trace the 
sources through many line-of-sight realizations. We find that, although the line-of-sight 
structure only occasionally changes the multiplicity of a family, the more common effect is to 
simply perturb the image positions. We find that these perturbations are also typically on the 
order of 1 arcsecond. 

We also study the effect of line of sight structure correlated with the cluster on the recovery of 
cosmological parameters using simulations. In our analysis, we utilize the multiple lensplane 
approximation. The total deflection is approximated as a series of point deflections occurring 
at discrete redshifts. The separation between lensplanes is typically ^100 Mpc. In other 
words, all matter within a rectangle of length ^100 Mpc is projected onto a plane. This is 
generally known to be a satisfactory approximation to treating the matter distribution as 
continuous. Regarding the effect of correlated large-scale structure, we note it can be 
reasonably modeled as sitting at the redshift of the cluster. Hence, it is reasonable to treat any 
correlated structure within ^50 Mpc as part of the projected mass of the lens itself. As an 
important example, we explored the effect of a long filament aligned along the line-of-sight 
behind the cluster. At the present, we are not able to utilize N-body simulations for this 
exploration. Instead, we use a simple parametric approach. We place a primary lens, 
representing a massive cluster (with mass M = 1.710^5 M©), at z = 0.184. We then place 5 
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equivalent circular potentials, with centers at the origin, spaced 10 Mpc apart behind the 
cluster to represent a filament (this is equivalent to breaking up the filament into 5 lensplanes). 
We use PIEMD profiles for these potentials with core radii of of 0.5 Mpc and cut radii of 1.5 
Mpc. With these parameters, each contributing potential has a constant density inside 0.5 Mpc 
and drops off as r"^ for r < 1.5 Mpc. Outside of 1.5 Mpc, the density goes as r"^. The total 
mass within the simulated filament is 10^^ Mq. We also create a "projected" model in which 
all potentials are placed on the same lensplane at z = 0.184. 

We lensed 10 simulated source catalogs with both the projected and full models and compared 
the resulting images. On average, deviations between the full and projected models are around 
0.6 arcsecond; less than half of the typical combined error from cluster galaxy scatter and 
uncorrected LOS halos. So even in the worst-case scenario of a long filament with mass 
roughly equal to the cluster, with its axis aligned exactly along the line-of-sight, the error due 
to modeling it in the same lensplane as the cluster is subdominant to the effects we have 
considered in our analysis. 

Bayesian MCMC approach to exploiting strong lensing 
observations 

We start with a parametrization describing the lens, the starting point as ddescribed earlier is 
the best-fit model for Abell 1689 pubhshed by Limousin et al. (2008). Using LENSTOOL 
software we explore the parameter space around the best-fit model while simultaneously 
solving for the family ratio for pairs of images. This is done while eproducing the location of 
the observed multiple images within the supplied uncertainties. Earlier versions of the 
software (Kneib et al. 1993) were based on a downhill minimization. However, that technique 
was very sensitive to local minima in the likelihood distribution; as a result, the modeling of 
complex systems as seen in Abell 1689 became rapidly too inefficient. In order to tackle the 
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current observational data for Abell 1689, we implemented a new optimization method based 
on a Bayesian Markov Chain Monte Carlo (MCMC) approach. The initial assumed model 
defines the prior PDF from which after sampling the posterior PDF is constructed for a range 
of acceptable model parameters. We use the Bayesian evidence to rank the models that best 
reproduce systems of multiple images. The Bayesian evidence is used to rank and characterize 
models when scanning parameter space. More details on this procedure can be found in Julio 
et al. (2008). 

In terms of the actual chains used, we run 10 interlinked Markov chains at the same time to 
prevent any Markov chain from falling into a local minimum. The MCMC convergence to the 
posterior PDF is performed with a variant of the "thermodynamic integration" technique 
(Ruanaidh and Fitzgerald 1996) called selective annealing. Here "selective" stands for the 
following process. At each step, 10 new samples (one per Markov chain) are drawn randomly 
from the current posterior PDF (which corresponds to the prior PDF at the beginning). These 
samples are then weighted according to their likelihood and selected with a variant of the 
Metropolis-Hasting algorithm (Metropolis et al. 1953; Hastings 1970). Roughly, the samples 
with the worst likelihood are deleted and the ones with the best likelihood are duplicated so 
that we always have 10 Markov chains running at the same time. These chains are run in 
parallel for efficiency. This ensures that the choice of initial seed does not impact the 
calculation, it generates more samples yielding smoother contour plots for the cosmological 
parameters and w^. 
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